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LONG-TERM  GOALS 

The  long  term  goal  of  this  effort  is  to  develop  a  methodology  to  estimate  bathymetry  and  wave-energy 
dissipation  from  satellite-based  SAR  imagery. 

OBJECTIVES 

This  program  will  attempt  to  develop  a  capability  for  estimating  bathymetry  via  assimilation  of  SAR 
imagery  into  a  near-shore  spectral  wave  model.  A  capability  for  estimating  dissipation  of  wave  energy 
due  to  bottom  friction  will  also  be  developed. 

APPROACH 

The  procedure  being  developed  makes  use  of  a  variational  data  assimilation  capability  for  the  SWAN 
near-shore  spectral  wave  model  of  Booij  et  al.  (1999).  The  bathymetry  estimation  procedure  makes 
use  of  the  SWAN  wave-action  balance  model  and  a  nonlinear  model  for  the  SAR  image  spectrum, 
along  with  adjoint  forms  of  both  models,  to  determine  the  bathymetry  which  minimizes  the  error 
between  the  estimated  and  observed  SAR-image  spectrum.  A  similar  approach  will  be  used  to 
estimate  wave  energy  dissipation  due  to  bottom  friction. 

WORK  COMPLETED 

The  work  completed  during  the  six  months  since  the  inception  of  the  program  includes  development  of 
the  mathematical  framework  for  the  assimilation  procedure  and  improvements  in  the  convergence 
properties  of  the  basic  SAR  assimilation  procedure. 

RESULTS 

The  expression  for  the  predicted  SAR-image  spectrum  derived  by  Hasselmann  and  Hasselmann  (1991) 
and  Krogstad  (1992)  can  be  written  in  a  slightly  modified  form  as 
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Si  {k,  ,ky)  =  jj G(x,  y,  , 

where 

G(x,y,kJ  =  {p^^{x,y)  +  f{x,y,k^)f{-x-y-k^)}e-^^^'^^'‘’^^ 
f{x,y,kj  =  l  +  ik^  [p^  ix,y)-  p„  (0,0)] ,  C{x,y)  =  p^  (0,0)  -  p„  {x,y) 

pM,y)  =  Re  jj|7;(^„^,)f 

PM,y)  =  ReJJ  T^(k^,ky)TXk,,ky)S(k^,ky)e‘^'^'‘^'^^^dkJky 

p„ix,y)  =  Re  jj|2;(A:„^^)f 

k'r(k„ky)  =  km{k^,ky)  and  T^{k k y)  = -^(gk^sinO- ia^  cosO). 

Here  S{kx,  ky)  is  the  wave  height  speetrum,  x  and  y  are  the  azimuth  and  range  eoordinates  and  kx  and  ky 
are  the  eorresponding  wavenumbers,  cr is  the  wave  frequency,  g  is  the  gravitational  acceleration,  R  is 
the  range  distance,  V  is  the  SAR  platform  velocity,  and  6  is  the  incidence  angle.  The  function  mikx, 
ky)  is  the  radar  modulation  transfer  function  (mtf),  which  describes  the  variations  in  the  radar 
backscatter  due  to  surface  tilt  and  hydrodynamic  effects. 

The  above  SAR  model  along  with  the  SWAN  model  form  the  basis  for  the  assimilation  procedure  used 
in  this  study.  The  methodology  proceeds  by  inputting  the  error  in  the  predicted  SAR  spectrum  into  the 
adjoint  SAR  model  (Lyzenga  2002)  whose  output  is  input  to  the  adjoint  SWAN  model.  The  resulting 
adjoint  wave  spectrum  can  be  used  in  a  procedure  to  minimize  the  error  between  the  predicted  and 
observed  SAR  image  spectrum.  Previous  implementations  of  this  assimilation  procedure  suffered 
from  poor  convergence  properties.  These  difficulties  have  been  traced  to  inaccuracies  in  the 
interpolation  procedures  used  to  go  between  the  polar  spectrum  produced  by  the  SWAN  model,  and 
the  rectangular  SAR  spectrum.  These  errors  would  cause  the  procedure  to  sometimes  diverge,  due 
primarily  to  lack  of  energy  conservation  in  the  polar  to  rectangular  conversion.  The  improved 
algorithm  now  converges  reliably. 

The  desired  bathymetry  estimation  algorithm  is  presently  being  implemented.  Estimation  of 
bathymetry  via  an  assimilation  approach  requires  a  first  guess  for  the  bathymetry  and,  at  each  iteration 
of  the  assimilation  loop,  a  bathymetry  correction  is  made.  The  depth-dependent  terms  in  the  SWAN 
model  form  the  basis  for  the  gradient  of  the  error  in  the  estimated  SAR  spectrum,  J,  with  respect  to  the 
bathymetry 

dJ  c dA  SC,  dASC  dA  SC^  dASC.]^  r  d  SS  ^ 

dh  ^  ydx  Sh  dy  &i  da  Sh  dO  Sh  j  ^  a  Sh 

where  A  (x,  s)  is  the  action  spectral  density,  H  (x  ,s)  is  the  adjoint  action  spectral  density,  and  the 
integrals  are  over  spectral  space.  Here,  the  SCj Sh  terms  represent  the  first  variation  of  the  given 
propagation  velocity  with  respect  to  the  water  depth  A(x),  calculated  from  the  expressions  for  the 
propagation  velocities  given  above,  and  the  SSj  Sh  is  the  first  variation  of  the  source  term  with  respect 
the  water  depth.  Using  this  expression,  a  spatially-varying  bathymetry  correction  can  be  calculated 
from  a  combination  of  the  predicted  wave-action  spectrum  and  the  adjoint  wave  action  spectrum, 
obtained  from  the  SWAN  model  and  its  adjoint. 
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IMPACT/APPLICATIONS 


This  represents  a  unique  approach  to  bathymetry  estimation.  Other  approaches  rely  primarily  on  wave 
kinematics,  and  use  the  dispersion  relation  to  determine  the  water  depth  for  the  observed  wavelength 
and  phase  speed.  In  the  above-described  approach,  the  variations  with  water  depth  of  wave  energy 
flux,  both  in  spectral  and  physical  space,  will  provide  the  bathymetry  estimate.  Hence,  the  spatial 
variations  in  the  wave  energy,  including  the  effects  of  refraction,  will  be  used,  rather  than  wave 
kinematics.  This  approach  does  involve  some  risks,  since  it  depends  on  the  accuracy  with  which 
spectral  and  spatial  wave  energy  distribution  can  be  calculated  using  the  assimilation  process  and, 
therefore,  to  some  extent  on  the  accuracy  of  the  modeling  of  the  source  terms  in  the  SWAN  model. 
However,  the  returns  could  be  great  if  accurate  bathymetry  estimates  can  be  obtained  over  an  entire 
100  km  X  100  km  area  from  a  single  ERS  SAR  image,  or  a  small  number  of  images. 
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